WIS-5/OO/Apr-DPP 
MAD-NT/00-04 
hep-ph/0004049 
April 2000 



Status of the solution to the solar neutrino problem 
based on non-standard neutrino interactions 

O" 

O ■ S. Bergmann^ M. M. Guzzo^ P. C. de Holanda^ 

■ P. I. Krastev^ and H. Nunokawa^ 

u 

Oh' 
< 



^ Department of Particle Physics 
Weizmann Institute of Science, Rehovot 76100, Israel 



^ Instituto de Fisica Gleb Wataghin 

T 1 

^ ' Universidade Estadual de Campinas, UNICAMP 

ON '. 13083-970 Campinas SP, Brazil 

. Department of Physics 

O '. University of Wisconsin, Madison, WI 53706, USA 

O 
O 

:a 

p^' Abstract 
> 

^ . We analyze the current status of the solution to the solar neutrino problem 

based both on: a) non-standard flavor changing neutrino interactions (FCNI) 
and b) non- universal flavor diagonal neutrino interactions (FDNI). We find 
that FCNI and FDNI with matter in the sun as well as in the earth provide a 
good fit not only to the total rate measured by all solar neutrino experiments 
but also to the day-night and seasonal variations of the event rate, as well 
as the recoil electron energy spectrum measured by the SuperKamiokande 
collaboration. This solution does not require massive neutrinos and neutrino 
mixing in vacuum. Stringent experimental constraints on FCNI from bounds 
on lepton flavor violating decays and on FDNI from limits on lepton univer- 
sality violation rule out i^e transitions induced by New Physics as a 
solution to the solar neutrino problem. However, a solution involving Vf, — > Vr 
transitions is viable and could be tested independently by the upcoming B- 
factories if flavor violating tau decays would be observed at a rate close to the 
present upper bounds. 
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I. INTRODUCTION 



In recent years the accuracy with which the solar neutrino flux is being measured has 
been improved significantly Better statistics and calibration of the pioneering experi- 

ments, as well as the first next-generation experiment SuperKamiokande, measuring the solar 
neutrino spectrum and the event rate as a function of the zenith angle with unprecedented 
precision, have provided a lot of new information about the solar neutrino problem 0. On 
the theoretical side several substantial improvements have been made in the standard solar 
model (SSM) ||7HT0| which now includes diffusion of helium and heavy elements and updated 
low energy nuclear cross sections relevant to the solar neutrino production [11|. Furthermore, 
the SSM has received an important independent confirmation by the excellent agreement 
between its predicted sound speeds and recent helioseismological observations 

All five solar neutrino experiments observe a solar neutrino flux which is smaller 
than predicted by the SSMs. In order to understand this discrepancy it has been suggested 
that neutrinos are endowed with properties which are not present in the standard electroweak 
theory |T^. These new properties allow the electron neutrinos to be converted along their 
way from the center of the sun to the detectors on earth into different neutrino flavors, i.e. 
into muon, tau, or possibly sterile [13] neutrinos. The fact that the terrestrial experiments 
are less sensitive to these neutrino flavors explains the observed lower counting rates. The 
most plausible solution is that neutrinos are massive and there is mixing in the lepton 



sector. Then neutrino oscillations in vacuum [jl4| or matter ||15| , p!6| (the Mikheyev-Smirnov- 
Wolfenstein (MSW) effect) can explain the deflcit of observed neutrinos with respect to the 
predictions of the SSM ^7^^. 

In his seminal paper Wolfenstein [IS] observed that non-standard neutrino interactions 
(NSNI) with matter can also generate neutrino oscillations. In particular this mechanism 
could be relevant to solar neutrinos interacting with the dense solar matter along their path 
from the core of the sun to its surface pO|-p^. In this case the flavor changing neutrino in- 
teractions (FCNI) are responsible for the off-diagonal elements in the neutrino propagation 
matrix (similar to the Am^ sin^ 2^ term induced by vacuum mixing). For massless neutri- 
nos resonantly enhanced conversions can occur due to an interplay between the standard 
electroweak neutrino interactions and non-universal flavor diagonal neutrino interactions 
(FDNI) with matter 



While many extensions of the standard model allow for massive neutrinos, it is important 
to stress that also many New Physics models predict new neutrino interactions. The minimal 
super symmetric standard model without i?-parity has been evoked as an explicit model that 
could provide the FCNI and FDNI needed for this mechanism. Systematic studies of the 
data demonstrated that resonantly enhanced oscillations induced by FCNI and FDNI for 
massless neutrinos pT|j25|] , or FCNI in combination with massive neutrinos l^l],^ can solve 
the solar neutrino problem. 

In this paper we investigate the current status of the solution to the solar neutrino prob- 
lem based on NSNI, which is briefly reviewed in section |I[ In the flrst part of our study 
(section |I11| ) we present a comprehensive statistical analysis of this solution. Our analysis 
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comprises both the measured total rates of Homestake , GALLEX , SAGE and Su- 
perKamiokande and, for the first time in the context of NSNI, the full SuperKamiokande 
data set (corresponding to 825 effective days of operation) including the recoil electron spec- 
trum and the day-night asymmetry. We have not included in our analysis the seasonal 
variation but we will comment on this effect. For the solar input we take the solar neu- 
trino fluxes and their uncertainties as predicted in the standard solar model by Bahcall 
and Pinsonneault (hereafter BP98 SSM) [^. The BP98 SSM includes helium and heavy 
elements diffusion, as well as the new recommended value [11| for the low energy S-factor, 
5*17 = 19^2 eV b. We also study the dependence of the allowed parameter space on the high 
energy neutrino flux, by varying the flux normalization as a free parameter. 

In the second part of our study for the first time a systematic, model-independent inves- 
tigation of the phenomenological constraints on FCNI and new non-universal FDNI relevant 
for solar neutrinos is presented (section [TVI ). Our two main goals are: a) to find out whether 
NSNI can be sufficiently large to provide a viable solution to the solar neutrino problem, 
and b) to study various kinds of new interactions in order to single out those New Physics 
models that can provide such interactions. Since the typical energy scales relevant for solar 
neutrinos are lower than the weak interaction scale and therefore lower than any New Physics 
scale, it is sufficient to discuss the effective operators induced by heavy boson exchange that 
allow for non-standard neutrino scattering off quarks or electrons. These operators are re- 
lated by the SU{2)l symmetry of the standard electroweak theory to operators that induce 
anomalous contributions to leptonic decays. Since SU{2)l violation cannot be large for New 
Physics at or above the weak scale, one can use the upper bounds on lepton flavor violating 
decays or on lepton universality violation to put model-independent bounds on the relevant 
non-standard neutrino interactions. 

We find that non-standard neutrino interactions can provide a good fit to the solar 
neutrino data if there are rather large non-universal FDNI (of order 0.5 Gi?) and small 
FCNI (of order a few times 1Q~^Gf)- Our phenomenological analysis indicates that FCNI 
could only be large enough to provide transitions, while v^. —>■ transitions are 

not relevant for the solution of the solar neutrino problem, because of strong experimental 
constraints. Large FDNI can only be induced by an intermediate doublet of SU{2)l (a 
scalar or a vector boson) or by a neutral vector singlet. We conclude that the minimal 
supersymmetric model with broken i?-parity |28| is the favorite model for this scenario. 

In section |V| we discuss how to to confirm or exclude the solution to the solar neutrino 
problem based on non-standard neutrino interactions by future experiments. We argue that 
the magnitudes of FCNI parameters necessary for z/g conversion in the sun could be 

tested independently by the upcoming _B-factories. Finally, we discuss briefly the possibility 
of distinguishing this solution from the others by future solar and long-baseline neutrino 
oscillation experiments. 
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II. NEUTRINO FLAVOR CONVERSION INDUCED BY NON-STANDARD 

NEUTRINO INTERACTIONS 



Any model beyond the standard electroweak theory that gives rise to the processes 

i^ef ^ i^ef , (1) 

l^af ^ l^af , (2) 

where (here and below) f = u,d,e and £ = fj.,T and a = e, /i, r, is potentially relevant for 
neutrino oscillations in the sun, since these processes modify the effective mass of neutrinos 
propagating in dense matter. 

The evolution equations for massless neutrinos that interact with matter via the standard 
weak interactions and the non-standard interactions in (|l]) and (@) is given by [pO| , ^] : 




^G,(?:^'l 5!/<:' 1 1 1 . (3) 




where Ae(r) and Ai{r) are, respectively, the probability amplitudes to detect a and I'i 
at position r. For neutrinos that have been coherently produced as in the solar core at 
position ro, the equations in @ are subject to the boundary conditions Ae(ro) = 1 and 
^^(^o) = 0. While ly-exchange of u^. with the background electrons gives rise to the well 
known forward scattering amplitude \/2GFn(,(r), the FCNI in (|lD induce a flavor changing 
forward scattering amplitude \/2GF4^nf{r) and the non- universal FDNI are responsible for 
the flavor diagonal entry y/2GF€'lnf{r) in eq. (H). Here 

""^^^^ \2n„(r)+n,(r) f = d 

is the respective fermion number density at position r in terms of the proton [neutron] 
number density Upir) [rinir)] and 

5 = e4^-^ and e' = e'l^^ ^^^^ (5) 

Kj F ^ F 

describe, respectively, the relative strength of the FCNI in (p, and the new flavor diagonal, 
but non-universal interactions in (||). G^^^,^ = r) denotes the effective coupling of 
the four-fermion operator 

Oi^{-^Vf,)[jj). (6) 

that gives rise to such interactions. The Lorentz structure of Ol depends on the New 
Physics that induces this operator. Operators which involve only left-handed neutrinos 
(and which conserve total lepton number V) can be decomposed into & {y — A)®{y — A) 
and a {y — A)®{y ^ A) component. (Any single New Physics contribution that is induced 
by chiral interactions yields only one of these two components.) It is, however, important 
to note that only the vector part of the background fermion current affects the neutrino 
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propagation for an unpolarized medium at rest [T^j2^. Hence only the {V — A) ^ (V) part 
of Ol is relevant for neutrino oscillations in normal matter. One mechanism to induce such 
operators is due to the exchange of heavy bosons that appear in various extensions of the 
standard model. An alternative mechanism arises when extending the fermionic sector of 
the standard model and is due to Z-induced flavor- changing neutral currents (FCNCs). For 
a discussion of Z-induced FCNC effects on solar neutrinos, see Refs. pO| , |3T| . 

A resonance occurs when the diagonal entries of the evolution matrix in eq. (Q) coincide 
at some point r^es along the trajectory of the neutrino, leading to the resonance condition 



(7) 



An immediate consequence is that new FDNI for / = e alone cannot induce resonant 
neutrino flavor conversions. 

As we will see in section |^ only conversions are compatible with the existing 

phenomenological constraints on el and e'{ . We note that in the minimal supersymmetric 



standard model with broken i?-parity |28] the relevant parameters are given by 



^331 ^131 



and 



I ^331 1 



1^131 1 



(8) 



in terms of the trilinear couplings X[j^ and the bottom squark mass M-^. 

The neutrino evolution matrix in eq. (§) vanishes in vacuum and is negligibly small 
for the matter densities of the earth's atmosphere. Therefore the probability of finding an 
electron neutrino arriving at the detector during day time is easily obtained by evolving 
the equations in (j^) from the neutrino production point to the solar surface. Furthermore, 
typically there are many oscillations between the neutrino production and detection point 
and a resonance. Therefore the phase information before and after the resonance is usually 
lost after integration over the production and detection region and one may use classical 
survival probabilities. Then at day time we have ETI 



pday 

J/e— >Ve 



\Mrs)f 



2 + ^2 



Pc)cos20^ cos 20^, 



(9) 



where is the solar surface position and in the analytic expression in eq. we denote 
by 6*^^ and 6*^, respectively, the effective, matter- induced mixing at the neutrino production 
point and at the solar surface. In terms of the New Physics parameters £, e' and the fermion 
densities the effective mixing is given by pO|j21[ 



tan2^„ 



(10) 



Note that tan26'm = 2e^ /(^'^ ~ 1) is constant for / = e. Pc is the level crossing probability. 



The approximate Landau- Zener expression is |20,21 



Pc = exp [—717/2] with 
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Up 



_d_ 

dx 



(11) 
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When neutrinos arrive at the detector during the night, a modification of the survival prob- 
abihty has to be introduced since the non-standard neutrino interactions with the terrestrial 
matter may regenerate electron neutrinos that have been transformed in the sun. Assuming 
that the neutrinos reach the Earth as an incoherent mixture of the effective mass-eigenstates 



Ui and z/2 the survival probability during night-time can be written as p2 



pnight ^ i^.%.-sin^g^ + P2e(l-2P,t^ 



(12) 

COS 2^^ ^ ^ 

Here is the probability of a transition from the state z/2 to the flavor eigenstate along 
the neutrino path in the Earth. 

For our analysis we assume a step function profile for the Earth matter density, which 
has been shown to be a good approximation in other contexts (see e.g. Ref. for a 
recent analysis of matter effects for atmospheric neutrinos). Then the earth matter effects 
on the neutrino propagation correspond to a parametric resonance and can be calculated 
analytically p^ . 



P2e = sin^ 9'^ + Wl cos 261 + W^W^ sin 29^ , (13) 

where the parameters Wi and contain all the information of the Earth density and are 
defined in Ref. [^. (The only difference is that in our case also the off-diagonal element 
of the neutrino evolution matrix varies when the neutrino propagates through the earth 
matter.) 

It is this interaction with the terrestrial matter that can produce a day-night variation 
of the solar neutrino flux and, consequently, a seasonal modulation of the data. (Note that 
this seasonal variation is of a different nature than the one expected for vacuum oscillations 
from the change of the baseline due to the eccentricity of the earth's orbit around the sun.) 



III. ANALYSIS OF THE SOLAR NEUTRINO DATA 

In this section we present our analysis of the solution to the solar neutrino problem based 
on neutrino flavor conversions induced by NSNI in matter. Our main goal is to determine 
the values of e and e' that can explain the experimental observations without modifying the 
standard solar model predictions. 



A. Rates 

First we consider the data on the total event rate measured by the Chlorine (CI) exper- 
iment [m, the Gallium (Ga) detectors GALLEX and SAGE ^ and the water Cherenkov 
experiment SuperKamiokande (SK) [^]. We compute the allowed regions in parameter space 
according to the BP98 SSM and compare the results with the regions obtained for an 
arbitrary normalization of the high energy neutrino neutrino fluxes. 

We use the minimal statistical treatment of the data following the analyses of 



Refs. |35,36]. Our x function is defined as follows: 
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ij=l,...,4 

where Rj^ and R"^'^ denote, respectively, the predicted and the measured value for the event 
rates of the four solar experiments {i = CI, GALLEX, SAGE, SK). The error matrix 
contains both the experimental (systematic and statistical) and the theoretical errors. 

In Fig. D the allowed regions in the parameter space of ef, and e'^ for neutrino scattering 
off c?-quarks are shown at 90, 95 and 99% confidence level (CL). In Fig. |l]a, the flux is 
fixed by the BP98 SSM prediction (/^ = 1). The best fit point of this analysis is found at 

et = 3.2 X 10"=^ and e't = 0.61 , (15) 

with Xmin — 2-44 for 4 — 2 = 2 degrees of freedom (DOF). Allowing an arbitrary flux 
normalization, a different best fit point is obtained for (ejj, e'^) = (2.2 x 10~^, 0.59) and 
/b = 1.36 with Xmin = 0-91 for 4 — 3 = 1 DOF. The result of this analysis is shown in 
Fig. |I|b. (Effects due to deviations of the hep neutrino flux from the standard solar model 
prediction are expected to be less significant and we do not consider them in this work.) 

In Fig. ^ the allowed regions in the parameter space of and e'" for neutrino scattering 
off M-quarks are shown at 90, 95 and 99 % CL. In Fig. the ®B flux is fixed by the BP98 
SSM prediction. The best fit point of this analysis is found at 

e!^ = 1.32 X 10"^ and e'l = 0.43 , (16) 



with Xmin — 2-64 for two DOF. Allowing an arbitrary ^B flux normalization Jb, a different 
best fit point is obtained for {e'^,eZ) = (5.8 x 10-^ 0.425) and Jb = 1.34 with xLn = 0-96 
for one DOF. The result of this analysis is shown in Fig. ^3. 

It is remarkable that the neutrino flavor conversion mechanism based on NSNI provides 
quite a good fit to the total rates despite the fact that the conversion probabilities (^ 
and ( p!2D do not depend on the neutrino energy. This is unlike the case of the vacuum 
and the MSW conversion mechanisms which provide the appropriate energy dependence 
to yield a good fit. For NSNI the only way to distinguish between neutrinos of different 
energies is via the position of the resonance r^es- Note that according to eq. (|^), Vres is a 
function of e' only. As can be seen in Fig. ^, n^/nf (/ = d,u) is a smooth and monotonic 
function of the distance from the solar center r, allowing to uniquely determine r^es for a 
given value of e' . From Fig. |^ it follows that a resonance can only occur if e'^ G [0.50, 0.77] 
for NSNI with d-quarks or G [0.40, 0.46] for NSNI with w-quarks. For both cases the 
major part of these intervals corresponds to r^es ~ 0.2 {Rq being the solar radius). For 
e'^„ within the 90% CL regions (indicated in Fig. |l] and Fig. |) we find Vres ~ 0.1 Rq. 
Since the nuclear reactions that produce neutrinos with higher energies in general take place 
closer to the solar center (see chapter 6 of Ref. [^] for the various spatial distributions of 
the neutrino production reactions), a resonance position close to the solar center implies 
that predominantly the high energy neutrinos are converted by a resonant transition. For 
r^es ^0.1 Rq practically all ^B- neutrinos cross the resonance layer, fewer ^Be- neutrinos pass 
through the resonance, while most of the pp-neutrinos are not be affected by the resonance 
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since their production region extends well beyond the resonance layer. Therefore for most 
of the allowed region in Fig. |I] and Fig. |^ we have 

(PCB)) < {PCBe)) < {Pipp)) . (17) 

We note that the above relation is still valid when taking into account that a significant 
fraction of the pp neutrinos crosses the resonance layer twice, if they are produced just 
outside resonance. This is - roughly speaking - because a u^, which undergoes a resonant 
flavor transition when entering the solar interior at r^gs is reconverted into a z/g at the 
second resonance when it emerges again from the solar core. In our numerical calculations 
we properly take into account the effects of such double resonances. 

An immediate consequence of the relation in eq. ( [T7| ) is that as long as /b = 1 the NSNI 
solution predicts that Rsk < Rci < Rca^ which is inconsistent with the observed hierarchy 
of the rates, Rci < Rsk < Rca, leading to a somewhat worse fit than the standard MSW 
solutions. However when treating fs as a free parameter, for fs ~ 1.3 — 1.4 the SK rate 
is sufficiently enhanced to give the correct relation between the rates. In this case also the 
neutral current contribution from z/^^^ e~ scattering is increased due to a larger u^^r flux, 
which is consistent with the Super-Kamiokande observations. We find that for the best fit 
points for {s,e') in Figs. and ^ and fs ~ 1-35 the survival probability for '^B, ^Be and 
pp-neutrinos are ~ 0.24, 0.4 and 0.7, respectively. 

In Fig. Pd and Fig. Qb the '^B neutrino flux normalization Jb has been varied as a free 
parameter in order to study the dependence of the allowed parameter space on the high 
energy neutrino flux. It is interesting to note that the allowed regions in these figures do 
not completely contain those in Fig. |I]a and Fig. ^a where the boron flux and its uncertainty 
is determined by the BP98 SSM. In order to explain this apparently inconsistent result we 
have have plotted ^ function of Jb in Fig. ^ allowing fs to vary within a sufficiently 

broad interval {0 < fs < 100) for every point in the {e, e') parameter space. The horizontal 
lines indicate the 68, 90 and 99 % CL limits for two DOF. The intersection of these lines with 
the curve determine the relevant ranges of the boron flux allowed by the experimental 
data. The vertical dotted lines indicate the la and 3a ranges of the boron neutrino flux in 
the BP98 SSM. 

Note that the minima are obtained for a boron flux significantly larger (/^ ~ 1.35) 
than the one predicted by the SSM (/^ = 1.0), as we already anticipated in the discussion 
of eq. (l7). Moreover, from Fig. ^ it follows that for < 1 the fit to the experimental data 



imposes stronger constraints on the boron flux than the SSM. For fs > 1 the situation is 
exactly the opposite. Therefore the effect of relaxing fs from its SSM value is that regions in 
the {e, e') parameter space where the averaged survival probability for ^B neutrinos, (P(^B)) 
is smaller can be easily compensated by a larger boron flux and obtain a lower value for xji- 
On the other hand regions where (P(^B)) is rather large require a small boron flux which 
is more difficult to achieve when eliminating the SSM constraint on Jb- This is the main 
mechanism behind the changes of the allowed regions upon relaxing the SSM constraint on 
Jb- It explains why the regions with large e are allowed in Fig. |l|a and Fig. ^a and are ruled 
out in Fig. |l|b and Fig. Qb. Here (P(^B)) is rather large and a small boron flux ~ 1 
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like in the SSM is preferred to explain the data. The opposite occurs in the area between 
the two disconnected regions in Fig. |I]a and Fig. Here (P(^B)) is comparatively small 
and therefore a larger boron flux increases Xr this region leading to the merging of the 
separated contours in Fig. |l|b and Fig. ||b when is treated as a free parameter. 



B. Zenith Angle Data 

Next, we consider the zenith angle dependence of the solar neutrino data of the Su- 
perKamiokande experiment. As mentioned above, NSNI with matter may affect the neu- 
trino propagation through the earth resulting in a difference between the event rates during 
day and night time. The data obtained by the SuperKamiokande collaboration are divided 
into five bins containing the events observed at night and one bin for the events collected 
during the day and have been averaged over the period of SuperKamiokande operation: 
403.2 effective days for the day events and 421.5 effective days for the night events. The 
experimental results suggest an asymmetry between the total data collected during the day 
(D) and the total data observed during the night (A^) 



N - D 

A = 2 = 0.065 ± 0.031{stat.) ± 0M3{syst.) . (18) 

In order to take into account the earth matter effect we define the following x^-function 
that characterizes the deviations of the six measured iZf"^) from the predicted (^*^) values 
of the rate function of zenith angle: 

Xl{e,e\az)= E ^ . (19) 

i=l,...,6 ^Z,i 

Here az,i refers to the total error associated with each zenith angle bin and we have neglected 
possible correlations between the systematic errors of these bins. Since we are only interested 
in the shape of the zenith angle distribution, we have introduced an overall normalization 
factor, azi which is treated as a free parameter and determined from the fit. (Using this 
procedure also prevents over-counting the data on the total event rate when combining all 
available data in section [HI E| ). Note that the experimental value of the day-night asymmetry 



in eq. (0) is not used in the fit, since the six zenith angle bins already include consistently 
all the available information about the earth effects. 

In Fig. ^ we show the allowed regions in the (e, e') parameter space for neutrino scattering 
off d- and u-quarks, respectively. The contours in Fig. ^ correspond to the allowed regions 
at 90, 95 and 99% CL. The best fit (indicated by the open circle) is obtained for (e^,e'^) = 
(0.251,0.62) and az = 0.819 with Xmm — ^-^O ^^r neutrino scattering off (i-quarks and 
at {tl, t'l) = (0.229, 0.690) and az = 0.685 with xLn = 1-44 for neutrino scattering off 
M-quarks (having 6 — 3 = 3 DOF in both cases). 

Finally, in Fig. ^ we show the expected zenith angle distributions for SuperKamiokande 
using the values of (e, e') determined by the best fit. For comparison, we also present in this 
figure the expected zenith angle distributions for the best fit values of (e, e') found in the 



combined analysis (that will be discussed in section [HI E|) . 
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C. Recoil Electron Spectrum 



We also consider the measurements of the recoil electron spectrum by Su- 



per Kamiokande The available data, after 825 days of operation, are divided into 18 

bins. 17 of these bins have a width of 0.5 MeV and are grouped into two bins for a super 
low energy analysis with energies between 5.5 MeV and 6.5 MeV and 15 bins with energies 
ranging from 6.5 MeV (the low energy limit) to 14 MeV. The last bin includes all the events 
with energies larger than 14 MeV. 

Since the electron neutrino survival probability does not depend on the neutrino energy 
in the NSNI scenario, the spectral distortion of the recoil electrons from neutrino due 
to the presence of a z/^,,- component in the neutrino flux is expected to be very small ||3^ 
and therefore, even a relatively small spectral distortion (such as the one expected in small 
mixing angle MSW solution) could rule out this solution. 

The x^-function that characterizes the deviations of the measured {S"^^) from the pre- 
dicted (5'*'^) values for the electron recoil spectrum therefore provides an important test of 
the NSNI solution. It is defined as: 



xl{e,e',as) 



i,j=l,...,18 



thi 









) 













thi 



obs 



(20) 



where the error matrix (squared) 



a. 



a1{stat.) + 



uncorr. , 



+ ai{corr.) aj{corr.) + aiitheor.) ajitheor.) (21) 



includes statistical [ai{stat.)] and systematic experimental errors (including both the uncor- 
related [ai{uncorr.y\ and the correlated [(Ti(corr.)] contributions) as well as the theoretical 



errors [aiitheor.)] (see Refs. [050 for more details.). Again, as in the analysis for the 
zenith angle dependence, we introduce an overall normalization factor a^, which is taken 
as a free parameter and determined from the fit, in order to avoid over-counting the data 
on the total event rate. Fitting the present data to our scenario we obtain Xmin — 20.0 for 
18 - 1 = 17 DOF, which is still acceptable at the 27% CL. 



D. Seasonal Variations 

The earth matter effects on neutrino flavor transitions induce a seasonal variation of the 
data (beyond the expected variation of the solar neutrino flux due to the eccentricity of the 
earth's orbit) due to the variation of the day and night time during the year. Since these 
variations can be relevant to other neutrino oscillation scenarios a positive signal could 



help to distinguish the various solutions and it is worthwhile to analyze the effects of such 
a variation in the NSNI scenario. 

The present SK solar neutrino data do not provide any conclusive evidence in favor of 
such a variation, but indicate only that the variation seems to be larger for recoil electron 
energies above 11.5 MeV. In our scenario, however, we do not expect any correlation between 
the seasonal variation and the recoil electron energies, since the electron neutrino survival 
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probability does not depend on the neutrino energy. Therefore any range of parameters 
that leads to a considerable seasonal modulation for energies above 11.5 MeV is disfavored 
by the data for lower energies. However, for the range of parameters {e, e') that can solve 
the solar neutrino problem, earth regeneration effects are never strong enough to induce a 
significant seasonal variation. Hence taking into account the data on seasonal variations 
neither changes the shape of the allowed region, nor the best fit points. 

E. Combined Analysis 

Our final result is the fit derived from the combined analysis of all presently available 
solar neutrino data. In Fig. |^ and Fig. ^ we show the allowed regions for (e^, e'^) and 
(e", e'^), respectively, using both the results from the total rates from the Chlorine, GALLEX, 
SAGE and Super Kamiokande solar neutrino experiments together with the 6 bins from the 
SuperKamiokande zenith angle data discussed previously. Although adding the spectral 
information to our analysis does not change the shape of allowed regions nor the best fit 
points, it is included in order to determine the quality of the global fit. However, we do 
not take into account the seasonal variation in our combined analysis, since the effect is 
negligible. 

For neutrino scattering off d-quarks the best fit for the combined data is obtained for 

4 = 0.028 and e't = 0.585 (22) 

with Xmin = 29.05 for 28 — 4 = 24 DOF, corresponding to a solution at the 22 % CL (see 
Fig. 0a). Allowing Jb ^ 1, the best fit is found at {et,e't) = (0.018,0.585) and fe = 1.38 
with Xmin = 26.62 for 28 — 5 = 23 DOF, corresponding to a solution at the 27 % CL (see 
Fig. ^). For neutrino scattering off w-quarks the best fit for the combined data is obtained 
for 

e"^ = 0.0083 and C = 0.425 (23) 

with Xmin = 28.45 for 28 — 4 = 24 DOF corresponding to a solution at the 24 % CL (see 
Fig. |a). Allowing Jb ^ 1, the best fit is obtained for (ej^, e'") = (0.0063, 0.426) and Jb = 1.34 
with Xmin = 26.59 for 28 — 5 = 23 DOF, corresponding to a solution at the 27 % CL (see 
Fig. |]b). These results have to be compared with the fit for standard model neutrinos, that 
do not oscillate (where the CL is smaller than 10"''), as well as to the standard solutions of 
the solar neutrino problem in terms of usual neutrino oscillations (36% CL) [|r^|TB|. 

Finally, in Fig. ^ we show the expected zenith angle distributions for SuperKamiokande 
using the best fitted values of (e, e') from the combined analysis. 

IV. PHENOMENOLOGICAL CONSTRAINTS ON e AND e' 

In this section we investigate whether the allowed regions for the parameters e^^ and e'l^ 
are at all phenomenologically viable. The analysis of non-standard neutrino interactions that 
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could be relevant for the solar neutrino problem is similar to the discussions in Refs. pO| , ^ 
where the possibility that FCNI explain the LSND results PSIBDI or the atmospheric neutrino 
anomaly |4^,^ was discussed. 



Generically, extensions of the Standard Model include additional fields that can induce 
new interactions: A heavy boson B that couples weakly to some fermion bilinears Bij with 
the trilinear couplings Xij, where i,j = 1,2,3 refer to fermion generations, induces the four 
fermion operator BjjBki at tree-level. The effective coupling is given by 



-iTl^' ^^^^ 

for energies well below the boson mass Mg. Thus, in terms of the trilinear coupling Aq,/ 
that describes the coupling of some heavy boson B to Ua (a = e, /i, r) and a charged fermion 
f = u,d,e the effective parameters in (13) are given by 



el = ^y:' and e'i = '^^^|! ~ . (25) 



Since any viable extension of the Standard Model has to contain the SM gauge symmetry, 
the effective theory approach presented in Refs. is completely sufficient to describe any 



New Physics effect for the energy scales typical to present neutrino oscillation experiments. 
Even though the effective theory obviously does not contain all the information inherent in 
the full high-energy theory, the parameters of the effective theory are all of what is accessible 
at low energies, when the "heavy degrees of freedom" are integrated out. 

The crucial point for our analysis is the following: Since the SM neutrinos are components 
of SU{2)l doublets, the same trilinear couplings Aq,/ that give rise to non-zero el^ or e'{^ 
also induce other four-fermion operators. These operators involve the SU{2)l partners of 
the neutrinos, i.e. the charged leptons, and can be used to constrain the relevant couplings. 

Noting that Lorentz invariance implies that any fermionic bilinear Bij can couple to either 
a scalar (5) or a vector (V) boson it is straightforward to write down all gauge invariant 
trilinear couplings between the bilinears (that contain SM fermions) and arbitrary bosons 
S and V that might appear in a generic extension of the Standard Model (see Tabs. 1-3 of 
Ref. [0). From these couplings one then obtains all the effective four fermion operators 
relevant to the solution to the solar neutrino problem in terms of NSNI as well as the 
S'[/(2)i-related operators that are used to constrain their effective couplings. (We do not 
consider here operators that violate total lepton number which can be induced if there is 
mixing between the intermediate bosons p^ .) 



While we refer the reader to Refs. PP],^ for the details of this model-independent 



approach, we present here two explicit examples relevant to solar neutrinos to demonstrate 
how SU{2)l related processes can be used to constrain the parameters el^ or e'l^. First, 
consider the bilinear Lfji (where L denotes the lepton doublet and f = e,u, d) that couples 
via a scalar doublet to its hermitian conjugate fnL. In terms of the component fields the 
effective interaction is 
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%^(^7^^/3) UriM - %^(^-7^^;3) UriM , (26) 



where = eL,fiL,TL for a = e,fi,T. A^/ is the trihnear couphng of LafR to the scalar 
doublet and Mi 2 denote the masses of its SU{2)l components. The important point is that 
the scalar doublet exchange not only gives rise to the four-Fermi operator Ol in (P) (with 
iy — A) ® iy + A) structure), but also produces the SU{2)l related operator 

0{^{Up)Uf), (27) 

which has the same Lorentz structure as Ol, with the neutrinos replaced by their charged 
lepton partners. Moreover, the effective coupling of Of , that we denote by G;^^, is related 
to Gl ,, by 

r2 



^i^y^ (28) 



2 



Constructing all the relevant four fermion operators that are induced by the couplings 
between the bilinears listed in Tabs. 1-3 of Ref . , one finds that in general Of is generated 
together with Ol' . Here /' can be different from / only for interactions with quarks, that 
is in some cases (Of) is generated together with 0^ (C?")- The leptonic operator Of is 
always generated together with Of, unless the interaction is mediated by an intermediate 
scalar SU{2)l singlet that couples to 

(L,L,), = -^(^ei-Ziz/,). (29) 

Note that the singlet only couples between two different flavors, since the coupling has 
to be antisymmetric in flavor space. Consequently a singlet that couples to the bilinear 
{LiLe)s cannot induce a non-zero ef^. The fact that the resulting four fermion operators 
only mediate FDNI is true because for the solar neutrinos we only care about v^. — > z/^ 
transitions. (For atmospheric neutrinos also — >■ u-r transitions induced by non-standard 
neutrino interactions with the electrons are of interest. In this case the coupling of {L^Le)s 
to (L^Le)J via singlet exchange inducing FCNI is possible p| .) 



The effective interactions that are mediated by a scalar singlet of mass M that couples 
to {LiLejs with the elementary coupling A^e are given by 



I A. 



M2 



I A. 



12 



M2 



where we used a Fierz transformation and the identity A'^'j^B'^ = —B'j^A to obtain (pOD. 
One can see that in this case Of is generated together with three more operators that have 
the same effective coupling (up to a sign). However, unlike for the case of intermediate 
doublets (or triplets), all these operators involve two charged leptons and two neutrinos. 
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A. Experimental constraints 



1. Flavor changing neutrino interactions 



There is no experimental evidence for any non- vanishing . Therefore, whenever (9/ 
is generated together with Ol, one can use the upper bounds on to derive constraints on 
Gl^^^ . The most stringent constraints on are due to the upper bounds on ji^ — > e~ e"*" e~ 
and T~ — > e~ e"*" e~ [|^,|^: 



BR(/i- 
BR(r- 



e e e 



< 1.0 X 10" 

< 2.9 X 10' 



-12 



(31) 
(32) 



Normahzing the above bounds to the measured rates of the related lepton flavor conserving 
decays, BR(/i~ e~ Ve ^fi) ~ 100 % and BR(r~ e~ Ue = 0.18 [||], we obtain 



< = G:,/Gf < 1.0 X 10- 
= G^^/Gi. < 4.2 X 10" 



(33) 
(34) 



Note that the bounds on e{ do only coincide with those for e^^ in the SU{2) symmetric 
limit. We will comment on possible relaxations due to SU{2)l breaking effects later in 
section [IV B . 

To constrain Gf we use the upper bounds on /i — e conversion from muon scattering 



off nuclei 46 



a{fi Ti — > e Ti) 



concluding that 



ct(/x- 


Ti - 


> capture) 




-Pb - 


e- Pb) 




Pb - 


> capture) 


a{fi 








32g _ 


-.Z/^32p*) 




4 


= Gl^/ Gf 



< 4.3 X 10 



< 4.6 X 10 



-12 



-11 



< 7 X 10 



-11 



(35) 
(36) 
(37) 

(38) 



is a conservative upper bound irrespective of the inherent hadronic uncertainties for such an 
estimate. 

To constrain G^^ we may use the upper bounds on various semi-hadronic tau decays that 
violate lepton flavor P5|,E6[: 



BR(r- 
BR(r- 
BR(r 
BRfr- ^ 



e-7r°) < 3.7 x 10"^ 

e"p°) < 2.0 X 10"^ 

■^e-f]) < 8.2 X 10"^ 

< 1.9 X 10"^ 



e TV^ TT 



(39) 
(40) 
(41) 
(42) 
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Let us first consider the tau decays into vr" and p°. Since these mesons belong to an isospin 
triplet we can use the isospin symmetry to normalize the above bounds ( |55| ) and (^) by the 
measured rates of related lepton fiavor conserving decays. Using BR(r~ i^tTt~) = 0.11 



and BR(r- ^ z/^p") = 0.22 we obtain 



Gl^in) < 8.2 X 10~^ Gf , and G^^(p) < 4.2 x 10"^ Gp ■ (43) 

Since the vr (p) is a pseudoscalar (vector) meson its decay probes the axial-vector (vector) 
part of the quark current. 

In general, any semi-hadronic operator Of can be decomposed into an / = and an 
1 = 1 isospin component. Only the effective coupling of the latter can be constrained by 
the upper bounds on the decays into final states with isovector mesons, like the vr and the p. 
If the resulting operator is dominated by the 1 = component, the bounds in (|43|) do not 
hold. But in this case we can use the upper bound on BR(r~ e~ rf) in (^). Since the rj is 
an isosinglet, isospin symmetry is of no use for the normalization. However, we can estimate 
the proper normalization using the relation between the r] and vr hadronic matrix elements. 



which is just the ratio of the respective decay constants, fr^/ f^c — 1-3 [^,|^. Taking into 
account the phase space effects, we obtain from ( p]) that 

6-^,(7/) < 1.1 X 10-2 Gf. (44) 

Since the 77 is a pseudoscalar meson its decay probes the axial-vector part of the / = 
component of the quark current, while the neutrino propagation is only affected by the 
vector part. As we have already mentioned, for any single chiral New Physics contribution 
the vector and axial- vector parts have the same magnitude and we can use ( ^4|) to constrain 



the isosinglet component of Of. In case there are several contributions, whose axial- vector 
parts cancel each other |5T], the / = component could still be constrained by the upper 



bound on BR(r- — > e~ ti^ n') in (^2]). While the calculation of the rate is uncertain due 
to our ignorance of the spectra and the decay constants of the isosinglet scalar resonances, 
we expect that the normalization will be similar to that of the vr, p and rj discussed before. 
Finally we note that the decay t" ^ e~ uo would be ideal to constrain the / = vector part, 
but at present no upper bound on its rate is available. 

While one can always fine-tune some parameters in order to avoid our bounds, our basic 
assumption is that this is not the case. Thus from (^31) and (HI) we conclude that 



el = GlJGF<lO-\ (45) 



2. Flavor diagonal neutrino interactions 

So far we have only discussed the upper bounds on FCNI. However, if the neutrinos 
are massless then in addition to the FCNI that induce an off-diagonal term in the effective 
neutrino mass matrix, also non-universal flavor diagonal interactions are needed to generate 
the required splitting between the diagonal terms. 
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In general any operator that induces such FDNI is related to other lepton flavor con- 
serving operators, that give additional contributions to SM allowed processes, and therefore 
violate the lepton universality of the SM. Then the upper bounds on lepton universality 
violation can be used to constrain these operators. Using the relation to the operators that 
induce the FDNI one may also constrain the latter. 

As we mentioned already for massless neutrinos only a non-zero e'l^ {q = u, d) can lead 
to a resonance effect, while FDNI that allow for scattering off electrons alone are insufficient 
to solve the solar neutrino problem. Therefore we only need to discuss the effective flavor 
diagonal operator 



(46) 



where a = e, fi,T and q = Ul^r, di^R. 

It is easy to check that for FDNI induced by heavy boson exchange Ol is always induced 
together with 



(47) 



where q' = UL,R,dL,R can be different from q. Moreover, intermediate scalar singlets and 
triplets (that couple to QL) as well as charged vector singlets and triplets (that couple to 
QL) also give rise to 



(48) 



where q' = u,d for q = d,u. (See ( |5DD as an example for FDNI mediated by an intermediate 
scalar singlet.) Since O^ induces an additional contribution to the SM weak decay tl ^ nur 
for a = r and to vr ^ Ua for a = e, fi, the relevant effective coupling can be constrained 
by the upper bounds on lepton universality violation in semi-hadronic decays. The latter 
leads to a deviation of the parameters 



1 r(7r- - 




^ 1 + 




A^r(7r- - 




Gp 


1 r(r- - 




^ 1 + 


G" - G'^ 


A^r(7r- - 




Gf 


1 T{W- 




^ 1 + 


Qq - 


NT{W- 


ez/e) 


Gp 



(49) 
(50) 
(51) 



from unity. Here denotes a normalization factor, which is just the ratio of the above two 
rates in the SM such that Ra/p = 1 if Gf^,^^ = Gl^^^. In the approximation we assume that 
Gl ^ Gp. From the most recent experimental data it follows that 



1.0017 ±0.0015 
1.005 ±0.005, 
0.987 ±0.023, 



(52) 
(53) 
(54) 
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implying that 

r'f _ r"i 

e7 = < 10^2 (55) 

Gr p 

is a conservative upper bound. If Of^, is induced together with O^, then in the SU{2)l 
symmetric hmit e'\ = e'l, but a modest relaxation due to SU{2)l breaking effects is possible 
(see section |1V B| ). 



It is essential to realize that not all New Physics operators that induce the FDNI relevant 
to solar neutrinos are related to Cf^,. For an intermediate SU{2)i scalar doublet (see eq. (^Gf ) 
for / = g) or a vector doublet (that couples to q^L) or a neutral vector singlet, only Of is 
induced together with O^. In this case one may only use the upper bound on Gfi that are due 
to the constraints on compositeness. The present data from pp —>■ e^e^ , fi'^fi" + X P9|j46 



imply an upper limit on the scale of compositeness A{qqlala) ~ 1.6 TeV, which translates 
into 

Gl, < lO-'Gp (56) 



as a conservative estimate for a = e, fi. (One-loop contributions to the Z width due to 0[ 



lead to a similar constraint on Gl^ IH^I-) However, no upper bound on A{qqTT) is available. 



For a neutral vector singlet Gf^i^ = Gl^^^ and from (|) and (|56|) it follows that 

e'X ^ 10^^ (57) 

while there is no model-independent bound on . For intermediate SU (2) l doublets 
the bound in (^Tj) could be relaxed somewhat, since the effective couplings of the relevant 
operators may differ due to SU{2)l breaking effects, which we discuss next. 



B. Constraining SU{2)l breaking effects 

The excellent agreement between the SM predictions and the electroweak precision data 
implies that SU{2)l breaking effects cannot be large. To show that the upper bounds on G^ 
(or Gi^) translate into similar bounds for if their related operators stem from the same 
SU{2)l invariant coupling, we recall from eq. (^) that in general the ratio of the couplings, 
^laUfi/^af3 (o^ GtaVfJ^Luis)-! is given by ratio M^/M|. Here Mi and M2 are the masses of 
the particles belonging to the SU{2)l multiplet that mediate the processes described by 
i^aup) respectively. If Mi 7^ M2 this multiplet will contribute to the oblique 



parameters |^T| S, U and, most importantly, T . A fit to the most recent precision data 
performed in Ref. determined the maximally allowed ratio (Mi/M2)^ax be at most 
6.8 (at 90 % CL) for intermediate scalars. (Vector bosons in general are expected to have 
even stronger bounds for the mass ratio). Consequently the upper limits on the effective 
couplings G^^ agree with those we derived for the corresponding G^ (or G^^) within an order 
of magnitude even for maximal SU{2)l breaking. Thus, barring fine-tuned cancellations. 

e4<6.8e{ and e'i^<6.8e'i, (58) 

at 90% CL. 
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V. IMPLICATIONS FOR FUTURE EXPERIMENTS 



In this section we discuss how to test the solution to the solar neutrino problem based 
on non-standard neutrino interactions in future neutrino experiments. 

Let us consider first the possibility of obtaining stronger constraints on New Physics 
from future laboratory experiments. Our phenomenological analysis shows that FCNI could 
only be large enough to provide z/g —>■ transitions while, model-independently, z/g — > 
transitions are irrelevant for solar neutrinos. Even for z/g transitions the required 

effective coupling has to be close to its current upper bound, which we derived from limits 
on anomalous tau decays. Therefore the solution to the solar neutrino problem studied in 
this paper could be tested by the upcoming i?-factories that are expected to improve the 
present experimental bounds on several rare r decays. For example, assuming an integrated 
luminosity of 30 fb~^ (corresponding to 3 x 10''r pairs) for the BaBar experiment, the 
upper limits on the branching ratios in (|39|)-(^2D could be reduced by one order of magnitude. 
This would decrease the bound on e% in ( ^Sj) to a value close to the smallest possible best 
fit values for e (c.f. Fig. |^ and Fig. |^) ruling out a large region of the parameter space and 
making the NSNI solution increasingly fine-tuned. 

Next we consider the implications for future solar neutrino experiments |53-57|. In Tab. 1 
we present the expected ranges for the event rates (normalized to the SSM expectation in 
the absence of neutrino fiavor transitions) of those experiments, if the solar neutrino problem 
is explained by NSNI. In Fig. ^ the predicted rates are presented graphically. The ranges 
correspond to the 95 % CL regions for [e, e') in Figs. ^ and ||b. As before we use the 
BP98 SSM predictions for the initial neutrino fiuxes and the survival probability in eq. (^ 
to compute the expected rates for each of the five detectors. Specifically, there are three 
types of detectors: (a) The Sudbury Neutrino Observatory (SNO) |5^, which is measuring 



the ^B neutrino charged current (CC) rate, (b) the BOREXINO |5|] and KamLAND |55 
experiments that are designed to observe the ^Be neutrino signal and (c) the HELLAZ ||56 
and HERON experiments dedicated to a precise measurement of the low-energy 
neutrino flux. 



Tab. 1: Future solar neutrino experiments and their rates predicted by the NSNI solution 



Experiment 


Start of operation 


Main neutrino source 


Rate predicted by NSNI 


SNO 


1999 


^B 


0.22 - 0.43 


BOREXINO 


2001 


^Be 


0.30-0.52 


KamLAND 


2001 


^Be 


0.30-0.52 


HELLAZ 


> 2002 




0.52-0.83 


HERON 


> 2000 


pp 


0.52-0.82 



The predictions for the rates in Tab. 1 reflect the relation between the predominant 
neutrino fluxes that we presented in eq. (p!?!), i.e. neutrinos with higher energies are in 
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general produced closer to the solar center and therefore more likely to pass through a 
resonance and undergo flavor conversion. 

As can be seen from Fig. ^ the suppression pattern of the NSNI solution is clearly different 
from the one predicted by the small angle MSW solution (c.f. Fig. 1 of Ref. [^). But there 
is a striking similarity between the NSNI solution and the LMA solution, including the 
preference for large fs (c.f. Fig. 7 of Ref. the absence of a spectral distortion and 

the modest day night effect. Consequently, using solar neutrino data, it will be difficult 
to distinguish the NSNI scenario from the LMA MSW solution. We note, however, that 
the KamLAND experiment will provide an independent test of the oscillation parameters of 
the LMA MSW solution by observing anti electron neutrinos from several nuclear reactors 
around the Kamioka mine in Japan. Thus, if KamLAND would indeed confirm the LMA 
MSW solution, then the NSNI solution discussed in this paper will be irrelevant. 



Since SNO already started taking data and is expected to have some results soon, 
let us consider some implication for this experiment. As we have pointed out one of the 
important features of the NSNI conversion mechanism is the absence of any distortion in 
the solar neutrino spectrum even though the averaged survival probabilities of neutrinos 
from different nuclear reactions in the sun are not equal. Due to this feature, the following 
simple relation between the SuperKamiokande solar neutrino event rate Rsk and the SNO 
CC event rate Rs^o (both normalized by the SSM predictions) holds: 

RsK = R^^oi'^-r) + rfB, (59) 

where Rsk and Rsno defined exactly as in eqs. (4) and (6) of Ref. [^] and r is given by 



J dEeR{E,) ^ 


f dE,<p'''{E,) 


^u^,re{Eu, Ee) 


j dE,R{E,) 


J dE,cp'^{E, 


)Cuee{Eu, Ee) 



^ = —r ^ \ ■ (60) 

Here, Ee and E^, are the electron and neutrino energy, respectively, R{Ee) is the Su- 
perKamiokande resolution and efficiency function, ^ is the ^B neutrino flux, and ai,^e and 
(Ju^^e denote the elastic scattering cross sections for z/gC" — ^ i^eC" and v^^t^-' — z/^,^e~, 
respectively. 

We note that Rsk and Rsno defined such that Rsk = Rsno the absence of 
neutrino flavor transitions including the case where 7^ 1. (Strictly speaking, a slight 
violation of the equality in eq. (|59D could be induced by the earth matter effect on these two 
experiments, since they are located at somewhat different latitudes.) Using the relation (p9D, 
the true flux of the ®B neutrino flux (0 ^)true = /s(0 ^)ssm could be precisely determined 
by combining SuperKamiokande and SNO solar neutrino measurements, if the solar neutrino 
problem is indeed due to NSNI. 

Finally let us discuss shortly the possibility of testing the solution studied in this paper 
by future long-baseline neutrino oscillation experiments. Since only Ue —>■ Vt transitions are 
viable, an independent test would require a {pT-) appearance experiment using an intense 
beam of z/g (z/g), which could be created at future neutrino factories (see, e.g., Ref. pOl). 
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Assuming a constant density and using the approximation that na — Uu — Sue in the 
earth, the conversion probabihty for a neutrino which travels a distance L in the earth is 
given by: 



36^2 + (1 -3e')^ 



sm 



36^2 + (1 - ^e'^V^GpneL 



(61) 



Numerically, the oscillation length in the earth matter can be estimated to be 



Lose ~ 8.1 X 10^ 



2 mol/cc 



1 



36e2 + (l-3£ 



7^2 



km. 



(62) 



Using eqs. (|6lD and (|62D and the approximation ng ~ 2 mol/cc (which is valid close to the 
earth surface), we find that, for the case of non-standard neutrino scattering off d-quark, 
P ~ few xlO"^ for K2K (L = 250 km) and P ~ few xlO'^ for MINOS (L = 732 km) for 
our best fit parameters. Similarly for w-quark, P ~ few xlO~^ for K2K and P ~ few xlO~^ 
for MINOS for the best fit parameters. These estimates imply that it would be hard but not 
impossible, at least for the case of scattering off (i-quarks, to obtain some signal of z/g — v^- 
conversions due to NSNI interactions by using an intense v^. beam which can be created by 
a muon storage ring [BO . 



VI. CONCLUSIONS 

According to our analysis non-standard neutrino interactions (NSNI) can provide a 
good fit to the solar neutrino data provided that there are rather large non-universal FDNI 
(of order 0.5 G^) and small FCNI (of order lO'^ - 10"^^^). The fit to the observed total 
rate, day-night asymmetry, seasonal variation and spectrum distortion of the recoil electron 
spectrum is comparable in quality to the one for standard neutrino oscillations. 

From the model-independent analysis we learn that NSNI induced by the exchange of 
heavy bosons cannot provide large enough v^. —>■ transitions, while u^. — Vj. FCNI in 
principle could be sufficiently strong. However, the current bounds will be improved by the 
up-coming i?-factories, providing an independent test of the NSNI solution. The required 
large non-universal FDNI (for z/g transitions into both and v^-) can be ruled out by the 
upper bounds on lepton universality, unless they are induced by an intermediate doublet of 
SU (2) L (a scalar or a vector boson) or by a neutral vector singlet. For — >• z/^ there exists a 
bound due to the limit on compositeness in this case, but for z/g — > z/^ there is no significant 
constraint at present. 

Generically only very few models can fulfill the requirements needed for the solution dis- 
cussed in this paper: massless neutrinos, small FCNI and relatively large non-universal 
FDNI. As for the vector bosons the most attractive scenario is to evoke an additional 
U^I^b-zLt gauge symmetry (where B is the baryon number and denotes the tau lep- 
ton number), which would introduce an additional vector singlet that only couples to the 
third generation leptons and quarks [0. Among the attractive theories beyond the stan- 
dard model where neutrinos are naturally massless as a result of a protecting symmetry, are 
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supersymmetric SU{5) models ||62| that conserve B — L, and theories with an extended gauge 
structure such as SU {3)c <S> SU (3) l ® U{1) n models where a chiral symmetry prevents 
the neutrino from getting a mass. These particular models, however, do not contribute 
significantly to the specific interactions we are interested in this paper. SU{5) models have 
negligible NSNI since they are mediated by vector bosons which have masses at the GUT 
scale. SU{3)c ® SU{3)l ® U{1)n models can provide large ee and e'g, but these models do 
not induce NSNI with quarks. From eq. (|^) it follows that no resonant conversion can occur 
in this case. 

Therefore we conclude that the best candidate for the scenario we studied are super- 
symmetric models with broken i?-parity, where the relevant NSNI are mediated by a scalar 
doublet, namely the "left-handed" bottom squark. Although in this model neutrino masses 
are not naturally protected from acquiring a mass, one may either evoke an additional sym- 
metry or assume that non-zero neutrino masses are not in a range that would spoil the 
solution in terms of the non-standard neutrino oscillations we have studied in this paper. 

Even though we consider the conventional oscillation mechanisms as the most plausi- 
ble solutions to the solar neutrino problem, it is important to realize that in general New 
Physics in the neutrino sector include neutrino masses and mixing, as well as new neutrino 
interactions. While it is difficult to explain the atmospheric neutrino problem and the 
LSND anomalies by NSNI 0,^], we have shown in this paper that a solution of the 
solar neutrino problem in terms of NSNI is still viable. The ultimate goal is of course a direct 
experimental test of this solution. The upcoming solar neutrino experiments will provide 
a lot of new information which hopefully will reveal the true nature of the solar neutrino 
problem. 



ACKNOWLEDGMENTS 



This work was partially supported by Fundagao de Amparo a Pesquisa do Est ado de Sao 
Paulo (FAPESP), Conselho Nacional de Desenvolvimento Cientffico e Tecnologogico (CNPq) 
and NSF grant PHY-9605140. We thank Y. Grossman, Y. Nir and C. Pena-Garay for helpful 
discussions. One of us (MMG) would like to thank the Physics Department, University of 
Wisconsin at Madison, where part of this work was realized, for the hospitality. 



21 



REFERENCES 



[1] T.B. Cleveland et al. (Homestake Collaboration), Astrophys. J. 496, 505 (1998). 

[2] W. Hampel et al. (GALLEX Collaboration), Phys. Lett. B 447, 127 (1999). 

[3] J.N. Abdurashitov et al. (SAGE Collaboration), Phys. Rev. C 60, 055801 (1999). 

[4] Y. Fukuda et al. (Kamiokande Collaboration), Phys. Rev. Lett. 77, 1683 (1996). 

[5] Y. Fukuda et al. (SuperKamiokande Collaboration), Phys. Rev. Lett. 82, 1810 (1999). 

[6] J.N. Bahcall, Neutrino Astrophysics, Cambridge University Press, 1989. 

[7] J.N. Bahcall and M.H. Pinsonneauh, Rev. Mod. Phys. 67, 781 (1995). 

[8] J.N. Bahcall, M.H. Pinsonneault, S. Basu and J. Christensen-Dalsgaard, 
Phys. Rev. Lett. 78, 171 (1997). 

[9] J.N. Bahcall, S. Basu and M.H. Pinsonneault, Phys. Lett. B 433, 1 (1998); see also J.N. 
Bahcall's home page, [http : / /www . sns . ias . edu/~jnb| . 

[10] See also, J.N. Bahcall and R.K. Ulrich, Rev. Mod. Phys. 60, 297 (1988); 
J.N. Bahcall and M.H. Pinsonneauh, Rev. Mod. Phys. 64, 885 (1992); 
S. Turck-Chieze et a/., Astrophys. J. 335, 415 (1988); 
S. Turck-Chieze and L Lopes, Astrophys. J. 408, 347 (1993); 

V. Castellani, S. Degl'Innocenti and G. Fiorentini, Astron. Astrophys. 271, 601 (1993); 
V. Castellani et a/., Phys. Lett. B 324, 425 (1994). 

[11] E.G. Adelberger et al.. Rev. Mod. Phys. 70, 1265 (1998). 

[12] For a detailed list of the references on these solutions see: 

"Solar Neutrinos: The First Thirty Years", ed. by R. Davis Jr. et al., 
Frontiers in Physics, Vol. 92, Addison- Wesley, 1994. 

[13] See for example: 

D.O. Caldwell and R.N. Mohapatra, Phys. Rev. D48, 3259 (1993); 

J. T. Peltoniemi, D. Tommasini, and J.W.F. Valle, Phys. Lett. B 298, 383 (1993); 

J.T. Peltoniemi and J.W.F. Valle, Nucl. Phys. B 406, 409 (1993); 

G. Dvali and Y. Nir, JHEP 9810, 014 (1998) ||hep-ph/9810257| ; 

G. M. Fuller, J. R. Primack and Y.-Z. Qian, Phys. Rev. D52, 1288 (1995); 

J. J. Gomez-Cadenas and M. C. Gonzalez-Garcia, Zeit. fur Physik C71, 443 (1996); 

N. Okada and O. Yasuda, Int. J. Mod. Phys. A 12, 3669 (1997), and references therein. 

[14] V. N. Gribov and B. M. Pontecorvo, Phys. Lett. B 28, 493 (1969). 

[15] L. Wolfenstein, Phys. Rev. D17, 2369 (1978). 

[16] S.P. Mikheyev and A. Yu. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985); 
Nuovo Cimento C9, 17 (1986). 



22 



[ir; 

[18 

[19 

[20 
[21 
[22 
[23 
[24 
[25 

[26; 

[27; 



J.N. Bahcall, P.I. Krastev and A.Yu. Smirnov, 
Phys. Rev. D 58, 096016 (1998) ||hep-ph/98072T6| ; 
Phys. Rev. D 60, 093001 (1999) |lhep-ph/9905220|| . 

M.C Gonzalez- Garcia, P.C. de Holanda, C. Pena-Garay and J.W.F. Valle, 
[hep-ph/9906469| , Nucl. Phys. B , in press. 

G.L. Fogli, E. Lisi, D. Montanino and A. Palazzo, [liep-ph/9912231j . 
M.M. Guzzo, A. Masiero and S.T. Petcov, Phys. Lett. B 260, 154 (1991). 
V. Barger, R.J.N. Phillips and K. Whisnant, Phys. Rev. D 44, 1629 (1991). 
E. Roulet, Phys. Rev. D 44, 935 (1991). 

S. Degl'Innocenti and B. Ricci, Mod. Phys. Lett. A 8, 471 (1993). 
G.L. Fogli and E. Lisi, Astroparticle Phys. 2, 91 (1994). 

P.I. Krastev and J.N. Bahcall, "FCNC solutions to the solar neutrino problem'' 
|hep-ph/ 9703267 . 



S. Bergmann, Nucl. Phys. B 515, 363 (1998) [|tiep-ph/970739S 

Resonant neutrino conversion induced by non-orthogonal massless neutrinos was first 
discussed by J.W.F. Valle in Phys. Lett. B 199, 432 (1987). 

However this mechanism can not induce a large effect on the solar neutrinos due to the 
stringent constraints on the model parameters. See also Refs. [^ ,^. 

[28] C.S. Aulakh and N.R. Mohapatra, Phys. Lett. B 119, 136 (1983); 

F. Zwirner, Phys. Lett. B 132, 103 (1983); 

L.J. Hall and M. Suzuki, Nucl. Phys. B 231, 419 (1984); 
J. Ellis et al, Phys. Lett. B 150, 14 (1985); 

G. G. Ross and J.W.F. Valle, Phys. Lett. B 151, 375 (1985); 
R. Barbieri and A. Masiero, Phys. Lett. B 267, 679 (1986). 

[29] S. Bergmann, Y. Grossman and E. Nardi, Phys. Rev. D 60, 093008 (1999) 
||hep-ph/9903517 [. 

[30] P. Langacker and D. London, Phys. Rev. D 38, 886 (1988); iMd 38, 907 (1988); 

H. Nunokawa, Y.-Z. Qian, A. Rossi and J.W.F. Valle, Phys. Rev. D 54, 4356 (1996) 
||hep-ph/960530i|] , and references therein. 

[31] S. Bergmann and A. Kagan, Nucl. Phys. B 538, 368 (1999) [ |tiep-ph/9803305|] . 

[32] S.P. Mikheyev and A.Yu. Smirnov, Sov. Phys. Usp. 30, 759 (1987). 

[33] M. Freund and T. Ohlsson, |hep-ph/990950ll . 

[34] E.Kh. Akhmedov, Nucl. Phys. B 538, 25 (1999). 



23 



[35] G.L. Fogli, E. Lisi and D. Montanino, Phys. Rev. D49, 3226 (1994). 
[36] G.L. Fogli and E. Lisi, Astropart. Phys. 3, 185 (1995). 

[37] Y. Suzuki, "Solar Neutrinos", talk given at the Lepton Photon Conference, 1999. 

See Fig. 8.2 of Ref. [^ for the recoil electron spectra shape from *B neutrino due to 

Ve.e~ — > Ve.t~ and t'^,re~ i^fj,,Te~ scattering. 

[39] P.C.de Holanda, CPena-Garay, M.C.Gonzalez-Garcia and J.W.F.Valle, Phys. Rev. 
D60 (1999) 093010. 



[40] S. Bergmann and Y. Grossman, Phys. Rev. D 59, 093005 (1999) ||hep-ph/9809524 l. 



[41] S. Bergmann, Y. Grossman and D.M. Pierce, Phys. Rev. D 61, 53005 (2000) 
||hep-ph/990939q| . 

[42] C. Athanassopoulos et al. (LSND Collaboration), Phys. Rev. Lett. 77, 3082 (1996); 
Phys. Rev. Lett. 81, 1774 (1998). 

[43] For recent analysis, see e.g.: 

M.C. Gonzalez-Garcia, et al, Phys. Rev. D 58, 033004 (1998); 

M.C. Gonzalez-Garcia, H. Nunokawa, O.L.G. Peres and J.W.F. Valle, 

Nucl. Phys. B543, 3 (1999) [|hep-ph/980730q] ; 



N. Fornengo, M.C. Gonzalez-Garcia and J.W.F. Valle, |hep-ph/0002147. 



[44] S. Bergmann, H.V. Klapdor-Kleingrothaus and H. Pas, |hep-ph/0004048 



[45] W. Bhss et al. (CLEO Collaboration), Phys. Rev. D 57, 5903 (1998) [|hep-ex/97120iq 



[46] Particle Data Group: C. Caso et al., Eur. Phys. J. C 3, 1 (1998); see also: 
[http://pdg.lbl.gov| . 

[47] J.F. Donoghue, E. Golowich and B.R. Holstein, Dynamics of the Standard Model, 
Cambridge University Press, 1992. 

[48] A. Pich, Review talk at Lepton Photon 99, Stanford University, August, 1999. 
See also: [hep-ph/9711279|, |hep-ph/9802257. 



[49] F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 79, 2198 (1997). 

[50] M.C. Gonzalez-Garcia, A. Gusso and S.F. Novaes, J. Phys. G, 2213 (1998) 
||hep-ph/9802254 ]. 



[51] M.E. Peskin and T. Takeuchi, Phys. Rev. Lett. 65, 964 (1990); 
Phys. Rev. D 46, 381 (1992). 

[52] The BaBar physics hook (SLAC-R-504), Chapter 12.2.6, available at: 

[http : //www, slac . Stanford, edu/pubs/slacreports/ slac-r-504.html| . 

[53] A. B. McDonald, for the SNO collaboration, Nucl. Phys. B 77 (Proc. Suppl.), 43 (1999). 



24 



[54] L. Oberauer, Nucl. Phys. (Proc. Suppl.) B 77, 48 (1999). 
[55] A. Suzuki, Nucl. Phys. (Proc. Suppl.) B 77, 171 (1999). 

[56] A. De Bellefon for HELLAZ collaboration, Nucl. Phys. (Proc. Suppl.) B 70, 386 (1999). 

[57] R. E. Lanou, in Proceedings of the 8th International Workshop on Neutrino Telescopes 
(Venice, Italy, 1999), ed. by M. Baldo Ceolin, Vol. I, page 139. 



[58] J.N. Bahcall, P.I. Krastev and A.Yu. Smirnov, [hep-ph/0002293 



[59] J.N. Bahcall, P.I. Krastev and A.Yu. Smirnov, Phys. Lett. B 477, 401 (2000) 
||hep-ph/99 1124811 . 

[60] S. Geer, Phys. Rev. D 57, 6989 (1998). 



[61] E. Ma and D. P. Roy, Phys. Rev. D 58, 095005 (1998) [|hep-ph/9806210 
E. Ma, Phys. Lett. B433, 74 (1998) ||hep-ph/970947l . 

[62] L.J. Hall, V.A. Kostelecky, S. Raby, Nucl. Phys. B B267 415 (1986); 
Y. Okada, |hep-ph/9809297| . 



[63] F. Pisano and V. Pleitez, Phys. Rev. D 46, 410 (1992); 
P. Frampton, Phys. Rev. Lett. 69, 2889 (1992). 



25 



FIGURES 



CO 



10" 



10-' 



10"' 



10"' 



Ga + CI + SK (Rates Only) 
(a) 'B flux fixed to BP98 SSM 



d-quark 




I I 90 % C.L. 

95 % C.L. 

99 % C.L. 



0.40 0.45 0.50 



0.55 0.60 
8' 



0.65 0.70 0.75 



10" 



10"' 



10"' 



10"' 



Ga + CI + SK (Rates Only) 



(b) B flux free 



I 90 % C.L. 
95 % C.L. 
99 % C.L. 



d-quark 




0.40 0.45 0.50 



0.55 0.60 

e' 



0.65 0.70 0.75 



FIG. 1. Region of e = ej^ and e' = e'jj which can explain the total rates measured by the Homes- 
take, GALLEX, SAGE and SuperKamiokande solar neutrino experiments in terms of non-standard 
neutrino interactions with d-quarks. (a) The best fit (indicated by the open circle) is obtained for 
= (0.0032,0.610) with xLn = 2-44 for 4 - 2 = 2 DOF. A second (local) ^ minimum (in- 
dicated by the solid square) is found at (e,e') = (0.034,0.610) with = 2.63. (b) Allowing for 
an arbitrary *B flux normalization /s, the best fit (indicated by the open circle) is obtained for 
(e,e') = (0.022,0.590) and /b = 1.36 with xLn = 0-91 for 4 - 3 = 1 DOF. 
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FIG. 2. Same as in Fig. |^ but for n-quarks. (a) The best fit (indicated by the open circle) is 
obtained for (e, e') = (0.0013, 0.430) with xLn = 2-75 for 4 - 2 = 2 DOF. A second (local) 
minimum (indicated by the sohd square) is found at (e, e') = (0.0083, 0.425) with = 2.70. (b) 
Allowing for an arbitrary fiux normalization fs, the best fit (indicated by the open circle) is 
obtained for (e,e') = (0.0058,0.425) and fs = 1-34 with xl^in = 0-96 for 4 - 3 = 1 DOF. 
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FIG. 3. The ratio of the number density of electron to that of d and u-quarks in the sun 
ne{r) /nf{r) (/ = d,u), is plotted as a function of the distance from the solar center. 
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FIG. 5. Region of e = ej^ and e' = e'^ which is excluded by day and night data (contained in 1+5 
bins) as measured by the SuperKamiokande solar neutrino experiment in terms of non-standard 
neutrino interactions with (a) d-quarks and (b) n-quarks. For d-quarks, the best fit (indicated 
by the open circle) is obtained for (e,e') = (0.251,0.620) and az = 0.819 with Xmin = ^-^O 
for 6 — 3 = 3 DOF. A second (local) minimum (indicated by the solid squared) is found at 
(e,e') = (0.0316,0.570) and az = 1.02 with = 5.20. For u-quarks, the best fit (indicated by the 
open circle) is obtained for {e,e') = (0.229,0.690) and az = 0.685 with xLn = 1-44 for 6 - 3 = 3 
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FIG. 6. Expected the zenith angle dependence with the our best fit values of (e, e') determined 
by the SK Zenith angle only as well as the combined analysis for (a) d-quarks and (b) u-quarks. 
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FIG. 7. The allowed region for e = ef, and e' = e'^ obtained by the combined analysis us- 
ing 4 rates + 6 zenith angle bins + 18 spectrum bins for non-standard neutrino interactions 
with d-quarks. (a) Fixing /b = 1 the best fit (indicated by the open circle) is obtained for 
= (0.028,0.585) with xLn = 29.05 for 28 - 4 = 24 DOF. There are two additional (lo- 
cal) minima at (e,e') = (0.0033,0.610) with = 29.40 (indicated by the solid square) and 
(e,e') = (0.21,0.61) with = 33.1 (indicated by the cross), (b) Same as in (a) but allowing a free 
Jb- The best fit (indicated by the open circle) is obtained for (e,e') = (0.018,0.585) and fs = 1.38 
with xLn = 26.62 for 28 - 5 = 23 DOF. 
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FIG. 8. Same as in Fig. but for u-quarks. (a) Fixing /_b = 1 the best fit (indicated by the 
open circle) is obtained for {e,e') = (0.0083,0.425) with xLn = 28.45 for 28 - 4 = 24 DOF. A 
second (local) minimum is found at (e, e') = (0.0013, 0.430) with = 30.27 (indicated by the 
solid square) (b) Same as in (a) but allowing a free Jb- The best fit (indicated by the open circle) 
is obtained for {e,e') = (0.0063,0.426) and /s = 1.34 with xliin = 26.59 for 28 - 5 = 23 DOF. 
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FIG. 9. The predicted ranges of the ratios of event rates for SNO, BOREXINO and HEL- 
LAZ/HERON to the corresponding event rates predicted from the SSM assuming that the neu- 
trino conversion induced by FCNI as well as FDNI is the solution to the solar neutrino problem. 
Histgrams indicate the predictions with the best fittes parameters whereas the error bars indicates 
the range determined by varying parameters within 95 % C.L. regions for (e,e') in Figs. |^b and 
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